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Abstract. Strongly coupled gravitational systems describe Einstein gravity and mat- 
ter in the limit that Newton's constant G is assumed to be very large. The nonlinear 
evolution of these systems may be solved analytically in the classical and semiclas- 
sical limits by employing a Green function analysis. Using functional methods in a 
Hamilton-Jacobi setting, one may compute the generating functional ('the phase of 
the wavefunctional') which satisfies both the energy constraint and the momentum 
constraint. Previous results are extended to encompass the imposition of an arbitrary 
initial hypersurface. A Lagrange multiplier in the generating functional restricts the 
initial fields, and also allows one to formulate the energy constraint on the initial hyper- 
surface. Classical evolution follows as a result of minimizing the generating functional 
with respect to the initial fields. Examples are given describing Einstein gravity in- 
teracting with either a dust field and/or a scalar field. Green functions are explicitly 
determined for (1) gravity, dust, a scalar field and a cosmological constant and (2) 
gravity and a scalar field interacting with an exponential potential. This formalism is 
useful in solving problems of cosmology and of gravitational collapse. 

PACS numbers: 0460, 9880H 



In the limit that Newton's constant G approaches infinity, the mathematical 
equations describing Einstein gravity coupled to matter simplify dramatically. One 
may safely drop second order spatial gradients while first order spatial gradients are 
retained. Using Hamilton-Jacobi (HJ) theory and the semidassical Green function 
method, a preceding paper (Salopek 1998) demonstrated how to construct a general class 
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of solutions for matter and gravity in the strong coupling limit. Here, those powerful 
methods will be extended to formulate the initial hypersurface problem. 

Strongly coupled gravitational systems are useful in studies of cosmology and of 
gravitational collapse. In cosmology, strongly coupled gravity has been used to describe 
the evolution of 'long-wavelength fields' arising, for example, from inflation which then 
serve as the primordial initial conditions for structure formation. The wavelength of 
these fields exceeds the local value of the Hubble radius, and different spatial points 
are no longer in causal contact (Salopek 1991, Salopek and Bond 1990 and Salopek 
1998). A strongly coupled expansion also appears in string theory formulations of 
cosmology (Veneziano 1997). In problems of gravitational collapse, strongly coupled 
gravity describes 'velocity-dominated' evolution where the gravitational potential terms 
which contain spatial gradients may be neglected as space collapses rapidly into a 
singularity. Numerical studies (Berger 1998, Berger and Moncrief 1993, Hern and 
Stewart 1998) indicate that many singularities that appear in general relativity are 
velocity-dominated. In addition. Hern and Stewart (1998) have shown numerically that 
the evolution of a certain class of Gowdy models may mimic that of velocity-dominated 
models even when spatial derivative terms are not negligible. 

The two essential ideas behind generating general solutions to strongly coupled 
gravity and matter are: (1) computing the semiclassical Green function solution of 
the energy constraint and (2) convolving the Green function with an arbitrary, gauge- 
invariant, initial state. In this paper. I will generalize the Hamilton- Jacobi formalism 
for strongly coupled gravity and matter to encompass the situation where one specifies 
the fields on some arbitrary initial hypersurface in superspace. 

Hamilton-Jacobi theory for general relativity is useful for two very important 
reasons: 

(1) It can be used to formulate a primitive theory of quantum gravity. At present, 
there are many problems in solving the Wheeler-DeWitt equation (DeWitt 1967) which 
is the canonical equation for quantum gravity. The worst is perhaps the problem of 
infinities. This issue may be temporarily postponed by expanding the wavefunctional 
in powers of h. At the lowest order, one must solve the Hamilton-Jacobi equation 
for general relativity which is mathematically self-consistent and does not require 
renormalization. This paper develops tools that allow one to solve the HJ equation 
in the strongly coupled limit. At higher order in h, infinities would certainly appear 
and they would have to treated at some later date presumably using some new ideas. 
However, an expansion in powers of h has been very successful in field theory, and one 
strongly suspects that it will be useful in the gravitational context although much work 
remains to be done. (Future prospects for a quantum description of cosmology and 
gravity have been discussed by Hartle 1997.) 

(2) It yields a covariant formulation of the gravitational field. Since its conception in 
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the 1960's, it has been known that the Wheeler- De Witt equation refers to the 3-metric 
7ab, but that it does not depend on the lapse N nor the shift N^. Quite often, this leads 
to problems in interpreting various approximate solutions. In a semiclassical setting, 
the Hamilton- Jacobi equation also does not depend on N and Ni. This is actually a 
blessing because a solution of the HJ equation is valid for all choices of the lapse and 
shift. In this sense, it provides a covariant description of gravity. 

HJ methods have been particularly useful in solving the long-wavelength problem 
of cosmology. In fact, they yield a transparent and elegant solution for the nonlinear 
evolution of long- wavelength fluctuations in cosmological models (Salopek 1991, Salopek 
and Bond 1990, Salopek 1998). (They are also being apphed to string cosmology 
by Saygih 1998). In the long- wavelength hmit, I feel that attempts to include 
nonlinear effects using higher order perturbation methods by Abramo et al (1997) 
are actually more difficult to apply and interpret, and they tend to obscure the 
fundamental simplicity of the long-wavelength problem: different spatial points evolve 
as homogeneous and independent universes. (Unruh 1998 has discussed the results of 
Abramo et al 1997 in more detail.) HJ methods are very powerful because they allow 
one to glue these independent spatial points using the momentum constraint to form one 
universe. This feature is crucial if one wishes to compute higher order terms in the strong 
couphng expansion (Parry, Salopek and Stewart 1994). Nonhnear long- wavelength fields 
also figure prominently in stochastic inflation (Vilenkin 1983, Starobinski 1986, Salopek 
and Bond 1991, Linde et al 1994, Vilenkin 1998). There, the probability distribution 
for long-wavelength fields is governed by a diffusion-type equation. 

Other applications of HJ theory to cosmology include: (1) the construction of 
inflation models which produce non-Gaussian fluctuations (Salopek 1992) — these 
models are still of observational interest (Fan and Bardeen 1992, Moscardini et al 
1993); (2) a relativistic formulation of the Zel'dovich approximation which describes 
how nonlinear pancake structures form in the distribution of galaxies (Croudace et al 
1994, Salopek et al 1994); (3) a covariant computation of density perturbations from 
inflation with one or two scalar fields (Salopek and Stewart 1995, Salopek 1995). 

In classic textbooks (Lanczos 1970, Goldstein 1981), one usually discusses Hamilton- 
Jacobi theory for unconstrained systems in the context of canonical transformations. 
However, for general relativity, the presence of constraints prevents a straightforward 
application of these standard methods. For example, one may have expected that one 
could apply the method of canonical transformations to gravity after a suitable phase 
space-reduction. However, since gravity is a nonlinear theory, this has not been achieved 
in practice except in some very special cases (linear perturbation theory, minisuperspace 
models, etc.). Rather, one should apply the powerful Dirac formulation of constrained 
systems to the constraints of general relativity (interacting with a scalar field just to 
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be specific), 

7-^[7r"^(x), 7r'^(x), 7afe(x), 0(x)] = 0, (energy constraint) (la) 

Hi[TT°'''{x),TT'^{x),^abix),(f){x)] = 0, (momentum constraint) (lb) 

where one replaces the momenta by functional derivatives of the generating functional, 
'S[yabix),(f){x)]: 

The two constraint equations ([T^b) are self-contained equations for S, and they may 
be taken as the starting point for a HJ formulation of gravity. If S is real, one recovers 
classical general relativity by integrating the definition of the momenta: 

{'lij - Ni\j - Nji,) /N = 2k 7-^/2 (2^,^^ _ ^. .^,;) __ , (3a) 

(0 - iV%) /AT = «: 7-1/2^, (36) 

where 

K = SnG = 8n/ml (4) 

denotes the gravitational coupling constant. Actually, HJ methods allows one to 
integrate the classical evolution equations (P^b) in an elegant way. 

In section 1, I give a brief review of the semiclassical Green function method of 
solving strongly coupled gravitational systems as advanced by Salopek (1998). In section 
2, I describe how to generalize the method by introducing a Lagrange multiplier in the 
generating functional to include the case of specifying the fields on an arbitrary initial 
hypersurface in superspace. In earlier work, it was found that one of the initial fields in 
the Green function must be set to zero. A specific choice was made previously, but here 
I would like to examine other possibilities. 

In section 3, I show how to reduce the number of degrees of freedom appearing in 
the energy constraint by four. In sections 4 and 5, I will construct nontrivial Green 
functions for the following situations: 

(1) Gravity, dust, a massless scalar field and a cosmological constant, 

(2) Gravity and a scalar field with exponential potential. 

Section 6 uses specific examples to illustrate the generalized Green function method 
as applied to gravity, dust and a scalar field (without cosmological constant). Section 
7 adds a cosmological constant. In section 8, I consider an advanced example involving 
gravity and dust. Here the initial hypersurface is defined in terms of the initial dust 
field and the Ricci scalar of the initial 3-metric. Since this system is not analytically 
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solvable, I expand the generating functional in a series. This example illustrates quite 
clearly how the energy constraint restricts the initial fields. I also show how to construct 
the 4-metric using the HJ formalism. A summary and conclusions are given in section 
9. 

(In this paper, the matter fields have been rescaled by factors of the gravitational 
coupling constant k, and consequently k disappears from most expressions. Consult 
Salopek (1998) for definitions.) 



1. Review of HJ solutions for strongly coupled gravity and matter 

In a HJ formulation, the energy constraint and the momentum constraint describing the 
3- metric ■yab{x) interacting with a dust field x{^) are, respectively, (see, e.g., Salopek 
1998), 

n{x)/K = 7"^/^ (27ac76d - tabled) T T h ^ = , (5a) 

7*fcT — + ^lk^ + ^x,i = o. (56) 
"^kjjj hki dx 

I will refer to {'^ab{x),x{^)) ^ the 'original fields'. The semiclassical Green function 
method of solving these two equations is described below. 



1.1. Semiclassical Green function method 



A complete solution of the energy constraint eq.(0) is one where the number of 
parameter fields equals the number of original fields. The Green function solution given 
below is a complete solution to the energy constraint with the additional property that 
the parameter fields, [7^°'' (x) , X(o) i^)] ) i^ay be interpreted as 'initial fields' on some initial 
hypersurface in superspace: 



{x{x) - X{Q){x] 



2y/S(o^)'cosh(W^z)-7 



1/2 



1/2 



where 



\^^r{\n[[h][h^^])\r.[[h][hl^)] 



(6 a) 



and [h] and are matrices, each with unit determinant, whose components are given 
by 



[h] 



ab 



(6c) 



6 



The integrand of the Green function, eg. (pa\j , has the famihar HJ form, 
1 

describing a free particle; the time t may be identified with the proper time x of the 
dust particle, and D"^ is the distance-squared between 'two points', {"Jab^llib), 
space of 3-metrics: 




This distance function, eq.([f^, was first derived by DeWitt (1967). 

In order to ensure consistency of the approach, it was necessary to set one of the 
parameter fields to zero. (This point will be discussed in further detail in section 2). 
For example, in Salopek (1998) the initial dust field was set to zero: 

X(o)(x) = 0. (8) 

One then constructs a general solution to the energy constraint and the momentum 
constraint through a superposition over an arbitrary 'initial state' X: 

S\lah{x),x{x)\ = G[lab{x),x{x)\ 1^! i^) , Xio){x) = 0] + X^i?] , (9a) 

where the initial 3-metric I'^abix) has been chosen to minimize S, 

Hb hab 

Here X is an arbitrary "gauge-invariant" functional of the initial 3-metric: 

That is, X[7^°''(a;)], is invariant under reparametrizations of the spatial coordinates. 
After solving the minimization equation (P^, one solves for the initial 3-metric as a 
function of the original fields, 

7i?(^)=7i?[7a.(x),x(x)] , (10) 

One then substitutes this result into eq.(^^ to determine the generating functional 
S[iab{,x),x{.^)] which depends solely on the original fields. 

The solution (|9^c) is an Ansatz. It was determined by trial and error. One can verify 
that it satisfies the energy constraint (|5^) by computing the functional derivatives of S 
with respect to jabix) and x{^)- K also satisfies the momentum constraint ( |5^ ) because 
the initial state X was chosen to satisfy the momentum constraint initially, eq. (|9^) , and 
the momentum constraint is preserved upon evolution because its Poisson bracket with 
the Hamiltonian vanishes. Many examples given in Salopek (1998) demonstrate how to 
exploit the Green function method and the superposition principle in practice. 



7 



1.2. Initial hypersurface prescription 

The simplest interpretation of setting the parameter field X(o){x) to zero, eq.(|^), is that 
one is choosing an 'initial hypersurface' where the original dust field vanishes. However, 
this choice is not forced upon us. In fact, in defining an initial hypersurface, one may 
set any one scalar function of the parameter fields to zero. I will refer to this process as 
the initial hypersurface prescription. For example, if, in addition to a dust field, there is 
a scalar field and if 0(o)(a;) denotes the corresponding initial field, one can instead 
set 

0(o)(x) = O. (11a) 

This choice corresponds to defining an initial hypersurface where the scalar field 
vanishes. In addition, one now assumes that X is a gauge-invariant functional, 

J = X[7l?(x),X(o)(x)], (lib) 

of the remaining parameter fields, 'y^^i^) X{o){x)- The initial hypersurface condition 
can presumably contain spatial derivatives such as, 

0?o)+7l?0Sy? = O, (lie) 

or even the Ricci scalar i?(o) associated with the initial 3-metric such as, 

X(o) — -B-R(o) = , (B is a constant). (Hc^) 

This last example will be treated in detail in section 8. 

In general one should be able to impose any one constraint on the parameter field 
such as eqs.( |lla| -c) by employing a Lagrange multiplier. I will develop this idea in the 
next section. 



2. The Lagrange multiplier method 

Although a Green function solution such as eq. (|6^ c) to the energy constraint (^^ is 
definitely a useful device, it leads to an inconsistency if one tries to construct solutions 
to the energy constraint and the momentum constraint by superimposing over all of the 
initial fields. This mathematical reason for this inconsistency is that the initial fields are 
not all independent since they obey the energy constraint on the initial hypersurface. 
Apparently, the initial energy constraint reduces the number of degrees of freedom by 
one field per spatial point. In order to avoid this problem, it was suggested in Salopek 
(1998) that one arbitrarily set one of the parameter fields to zero before superimposing 
over the remaining parameter fields. As was argued on physical grounds in section 1.2, 
this choice defines the initial hypersurface. In the present section, I will expand on this 
theme by introducing a Lagrange multiplier L{x) to impose this reduction in the number 
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of parameter fields. This simple extension leads to a much wider field of applications 
for the HJ formalism. The role of the initial energy constraint will be discussed using 
an example in section 8. 



2.1. Formal HJ solution for gravity, dust and a scalar field 

In this section, I will construct a formal HJ solution to a strongly coupled system 
consisting of gravity, a dust field and a scalar field (f){x) with potential V{(f)). This 
system will be described by the following constraint equations: 

, , , -1/2 /fA \ '^'5 

n{x)/K= 7 ' {2-faclbd - lablcd) T h ^ 

O'Jab O'Jcd OX 

n, , . f SS \ 6S SS ^ SS , 

7^i(x) = -2 7ifc- — + ^7fc/,* + T- 0,* + ^ X,* = . (125) 
V o-fkjjj d-fki d(p dx 

For this case, there are eight degrees of freedom per spatial point associated with the 
original fields, •jabix), </)(x) and x{x). 

An Ansatz solution to the energy and momentum constraints is 

S[^ab{x)Aix),xix)]=g + T + C. (13) 

I will assume that the Green function Q satisfies the energy constraint eq.( [L2c| ) with S 
replaced by Q: 

5G sg sg 



n{x)/K= 7 (27ac7M - 7a67cd) 



hab hcd 5x 



1 



-1/2 fsg^ 2 



, 2 yj^) + V^^n^) = 0. (14a) 

It is a complete solution in that it depends on eight inhomogeneous parameter fields 
(initial fields), Jabix), 0(o)(a;) and X{o){x)- 

g = g[lab{x),(j){x),x{x)\labi^)^<P{0)i^)^X{0){x)] ■ (146) 

I, the 'generalized initial state', is a functional of all the parameter fields, 

^ = AlS! i^)^ 0(0) (x) , X{o) (x)]. ( 14c) 

The new ingredient, C, is a functional of the parameter fields which is linear in the 
Lagrange multiplier L{x), 

^ = J d^x-fl/^^ L{x) f [-fabix), (l)(o)ix), X(o)ix)] . (Ud) 

Here / is some scalar function of the parameter fields, 7ife''(a;), (f)(Q){x) and Xio){x) which 
will, in general, contain spatial derivatives. I will demand that X and C be invariant 
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under spatial coordinate reparameterizations of the parameter fields and the Lagrange 
multiplier L{x): 



,(0) \ _^ JO) . AO) . (0) 



o = -2 hr-i) +Tf^0r + T^xr as^) 

^tSV,,- ^Tki ^^(0) ^Xio) 

^Tkjjj STm ^'^(o) '^^(0) 

The parameter fields and the Lagrange multiplier are chosen to minimize the 
generating functional, eq.(0): 

(50(0) (x) 

= -^^, (16c) 

The minimization process leads to a solution of the classical Einstein equations in the 
strongly coupled limit. 

Variation with respect to the Lagrange multiplier L{x\ eq. QlGtj) , implies that the 
parameter fields are constrained according to, 

= / [7l?(a;),0(o)(a;),X(o)(a;)] • (17) 

This effectively reduces the number of parameter fields by one in accordance with 
the initial hypersurface prescription of section 1.2. The minimization method will 
be successful if one can solve eqs.( |16a| -d) for the parameter fields and the Lagrange 
multiplier in terms of the original fields: 

itK^) = lab [lah{x), 0(X), X{X)] , (18a) 

0(0) (x) = 0(0) [7„b(x),0(x),x(a;)] , (18&) 
X(o)(a;) = X{o) [labix), 0(a;), xix)] , (18c) 
L{x) = L [7„b(x),0(x),x(x)] . {18d) 

If, for example, there are too many parameter fields, something will go wrong at this 
stage. Substitution of the above into eq.(|13|) leads to the ultimate goal, the generating 



functional S, which depends only on the original fields. 
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2.2. Justification of Ansatz 



In order to justify the Ansatz (|T3|), it is necessary to compute the functional derivatives 
of the generating functional with respect to the original fields. 

Before one applies the minimization prescription, iS is a functional of the original 
fields, the parameter fields and the Lagrange multiplier: 



After applying the minimization prescription, the parameter fields and the Lagrange 
multiplier become functions of the original fields through eqs.(^8^d). Hence a functional 
derivative of S with the respect to one of the original fields, say can be computed 
by the chain rule: 



6S 



6(l){x) 



5g 



7a6,X 



6(l){x) 
6S 



+ 



7a6,X,7i°\</'(0).X{0) 



6S 



+ 



5S 50(0) (y) 



+ 



6S 



50(0) (y) 50(a;) 
5Liy) 



+ 



(20 a) 



Sxmiy) Hix) SL{y) 6(l){x) 

All of the terms in eg. ( ^Ooj ) which are integrated over y vanish by virtue of the 
minimization conditions, eqs.( |16Q| -d). Hence a functional derivative of S with respect 
to 4>{x) coincides with a functional derivative of Q with respect to 0(x): 



6S 



50(a;) 



5g 



7a6,X 



6(l){x) 



This argument also holds for functional derivatives with respect to the other original 
fields, 'jab{x) and x(x). 

Hence, S given by eq.(|13D satisfies the energy constraint (|12q|) because: 

(1) Q satisfies the energy constraint ( |14aD , and 

(2) functional derivatives of Q and S are identical, which was illustrated in the previous 
paragraph. 

Moreover, S also satisfies the momentum constraint ( |12&|) since I demanded 
invariance under reparameterization of the spatial coordinates of both X and C through 
eqs.( p.5Gj -b). Since reparametrization invariance is valid initially, it is guaranteed to be 
valid at other times. 

Thus, S given by eq.(p!3D is a HJ solution to both the energy constraint and the 
momentum constraint provided one can solve for the initial fields and the Lagrange 
multiplier in terms of the original fields, eqs.(^8^d). 



(206) 
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2.2.1. Discussion of Lagrange multiplier method I will now demonstrate the physical 
significance of the Lagrange multiplier method. The functionals, X and C, appearing in 
eq.(|T3D effectively define the 'initial setting' of the parameter fields. 

As was mentioned briefly in section 1.2, different choices of the / function, eg. (p!7D , 
specify the initial hypersurface. For example, if 

/ = X(o), (21a) 

then variation of the Lagrange multiplier implies that X{o){^) vanishes. After setting 
X{o){x) = in the generating functional, eq.(|T3D, one recovers eq.(|9^c) which is the 
result given in Salopek (1998). There, it was shown by explicit construction that if the 
dust field were zero, x{^) = 0, then the generating functional S and the initial state I 
coincide: 

Sbabix), xix) = 0] = X[7l?(x) = labix)] . (215) 

In this sense, x{^) = denotes the initial hypersurface which was chosen among all 
possible choices by setting the parameter field X(o){^) = 0. 
If on the other hand, 

/ = 0(0) , (22a) 

then the scalar field is uniformly zero on the initial hypersurface: 0(x) = 0. As was 
shown in Salopek (1998) for the case of a single scalar field interacting with gravity, if 
4>{x) = 0, the generating functional S and initial state X coincide: 

S[^ab{x), 0(x) = 0] = liiSlix) = 7afe(x)] . (225) 

Once again, setting a parameter field to zero effectively determines the initial 
hypersurface in superspace. 

The choices for / are limitless and they really depend on the problem at hand. 
The reader should be warned, however, that if / depends on spatial derivatives of the 
parameter fields, then it may be difficult in practice to solve for the Lagrange multiplier 
and the parameter fields as required by eqs. (|18Q| -d). Such a case is discussed in section 
8. 

The general formalism is now complete. The basic result is encapsulated in the 
simple expression for the generating functional, eq. ([T3|) : S = Q + I + C The 
Green function Q effectively describes how the system described by the original fields 
['yab{x),4'{x),x{x)] evolves from initial fields [7^°''(x), 0*^°^(x), x^°''(a;)]- X and C describe 
the initial setting. X is the initial state; £, the Lagrange multiplier term, prescribes the 
initial hypersurface. Classical evolution follows from the minimization of the generating 
functional with respect to the initial fields and the Lagrange multiplier. 
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3. Reduced energy constraint 

In sections 4 and 5, I will expand the list of known Green functions to include (1) gravity, 
dust, a massless scalar field and a cosmological constant and (2) gravity and a scalar 
field with an exponential potential. Before constructing these solutions, I will describe 
how to reduce the energy constraint to a more manageable form. 

In searching for the Green function solution to the energy constraint eq. (|12Q| ) 



describing gravity, dust, and a scalar field with potential, one must deal with eight 
degrees of freedom: two matter fields and and the six components of the symmetric 
matrix, 'jab- One may think that it is hopeless to solve this equation, but the presence 
of symmetries makes the task manageable. One attempts a solution of the form, 

g = g[(f){x),xix),cyix),z{x)], (23a) 

where the 3-metric degrees of freedom are parameterized by a and z, 



a = ^ln7, z= yTT{\n{[h][hl^^])\n{[h][hi;^])}, (236) 



and [h] and [h^^^] were defined in eq.(|6^. The 'reduced energy constraint' involves only 
four degrees of freedom: 

(reduced energy constraint). (23c) 



Not every solution to the energy constraint has the form ( p3a|) but the Green function 



does. In fact, the Green function solutions given in sections 4 and 5 are derived from the 
reduced energy constraint (|23(j ). It is much easier to work with the reduced equation 
than the full energy constraint ( P.2g| ) because the former does not contain tensor indices. 

4. Green function solution: gravity, dust, massless scalar field with 
cosmological constant 

If the potential V{(j)) = Vq describes a cosmological constant, then the Green function 
solution of the energy constraint ( |12c| ) for gravity, dust and a massless scalar field is: 



f 2,H 

G[lab{x),x{x),(p{x)\labix)^X{0){x),(l){0){x)] = / d^X 



sinh 



3Ho 
2 



(x - X(0) 



|27^/Sf/;cosh [^f^^ z- + {<!>- 
- (7'^' + 7('o)') cosh ^ (x - X(o)) 



(24a) 
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where 




(245) 



This solution depends on eight parameter fields, (^7^^ (x), X(o)(a;), 0(o)(a;)j. I will now 
discuss special cases of this rich solution. 



4.1. Gravity, dust and cosmological constant 



The Green function solution for gravity, dust and a cosmological constant {i.e., no scalar 
field) is found by setting 



(25) 
in eq.( p4aD . A complete solution depending on seven parameters, [^^^^ {x) , X{o){^)) 1 is: 

2H, 



sinh 



2y/Sf^coshL/^.] -(7^/^ + 7?^) cosh ^(x-X(o)) 



(26) 



One can readily verify that the above solution reduces to the case of gravity and dust, 
eq.(^), in the limit that the cosmological constant vanishes, Hq 0. 

In general, it is not permissible to arbitrarily specify original fields in the Green 
function ( p4a[) , and then expect that they will satisfy the energy constraint. However, 
the solution given in eq.(^) may be derived from eq. ( p4aD using the following argument. 
If one assumes that the functionals, X and C, in eq.(0) that define the initial setting are 
independent of 0(o)(x), then eq.(^) follows from minimization with respect to 0(o)(a;). 
As a result, the scalar field dependence in the Green function ( P4q| ) drops out, and one 
loses this degree of freedom. 



4.2. Gravity, massless scalar field and cosmological constant 
The Green function, 



G[lab{.x),(t){x) \ 7!? (a;), 0(0) (a;)] 



/4H 



d'^x 



7 + 7(0) - 27^/27(0/)^ cosh W -y/^^^ (0-0(0)) 



1/2 



(27) 



for gravity, a massless scalar field and a cosmological constant {i.e., no dust field) was 
given in Salopek (1998). The sign preceding the Green function is arbitrary and I have 
chosen it to correspond to an expanding universe. One may derive it from the Green 
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function eg. ( |24q| ) in the following way: if the functionals, X and £, in eq.(0) that define 
the initial setting are independent of X{o){^)y then minimization of S with respect to 
X{o)ix) implies 

sg 



which yields eg. (|27|) . Thus the dust degree of freedom drops out. 



(28) 



4.3. Gravity, dust and a massless scalar field 

The Green function for dust, gravity and a massless scalar field [i.e., no cosmological 
constant) is found from eq.( ^4QD by considering the limit that the cosmological constant 
vanishes, Hq ^ 0: 



4 /" 1 

g[-fab{x),(j){x),x{x)\-fi!j!{x),(j)(Q^{x),X{o){x)] = - d^x ^— — 

'^•^ [X{x) - X(o){x, 



27^/^ 7(0^; cosh 



1/2 
^(0) 



(29) 



4-4- Special Case: 7(0) — >■ limit 

Note that in the limit that the determinant of the initial 3-metric approaches zero, 
7(0) 0, the Green function for gravity, dust, scalar field and a cosmological constant, 
eq.( |24Q| ), takes on the simple form 

J d'x^'/^ h{x\X{0)) , (30a) 
which describes the integral over the original 3-geometry of the the Hubble function 



H (xlX(o)) = i^ocotanh 



3H, 



- [X- X(o) 



(306) 



The Hubble function describes gravity, dust and a cosmological constant (the scalar 
field drops out), and it was first computed by Salopek and Stewart (1992). 

5. Green function solution: gravity and scalar field with exponential 
potential 



Consider a scalar field with exponential potential 

VP/ 

The reduced HJ equation (E^ becomes: 



V(0) = roexp 



-12^ 



-3a 



6a 



-3a 



5z 



' 1 



(31) 



3a 



+e=^"K)exp ( 
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The Green function solution which depends on seven parameter fields [7ifc''(a;), 0(o)(a;)] 
is 



^=-\/3(l-V(3p))/'''^V7- 
where / is an abbreviation for 

p "X l'^ I -\ l't>{0) „ 1/2 1/2 

/ = 7 e V +7(0) e y " -2-^ ' 7^0') exp 



and 



e 



2\ 



n 2 



ln(7/7(o)) 



+ 1 



3p 



;z2 . 



(33 a) 



cosh ^,(33 6) 



(33 c) 



The sign in front of the Green function is arbitrary and I have chosen a negative sign 
to describe an expanding universe. 



5.1. p 



00 



As p — > 00, the potential for the scalar field approaches a constant, and one recovers the 
case of gravity, a massless scalar field and a cosmological constant described by eg. (pTj) . 



5.2. Special Case: 7(0) —* limit 

Note that as 7(0) — ^ 0, the Green function approaches 

Q -^-'^ j rf'x7^/2 i/(0), 
where the Hubble function was computed by Salopek and Bond (1990): 

H{4>)- 



3(l-l/(3p)) 



exp 



2p^ 



(34 a) 



{Uh) 



6. Examples of Lagrange multiplier method: gravity, dust and scalar field 

1 will illustrate the Lagrange multiplier method by considering various cases of gravity, 
dust, and a massless scalar field whose evolution is described by the Green function 
eq.(|29|). The case with cosmological constant is also interesting but the algebra is more 
complicated and I will devote section 7 to it. 

The Lagrange multiplier method for solving the HJ equation had been foreshadowed 
in an earlier paper (Salopek 1991). I will show how to recover these earlier results from 
the generalized formalism described in section 2. 

In the present section, I will consider the special case where the functionals X and C 
which specify the initial setting in eq.(p!3|) do not contain any spatial gradients. Many 
examples are given. 
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6.1. Initial setting contains no spatial gradients 

I will choose the function / appearing in the Lagrange multiplier term ( |14dD to be an 
arbitrary function of 0(o) and X(o)- 



/ = /(0(o),X(o)^ 
and I will choose the initial state 



J 



(35) 
(36) 



to be an integral over the initial volume of some arbitrary function g of 0(o) and X(o)- 
This system will be tractable because no spatial gradients appear. 

Minimization of the generating functional eq. (|T3|) with respect to the initial 3- metric 
labi^) implies 



1/4 
^(0) 



0, 

7^^^ cosh 



4 



{x - xm) {g + Lf) 



(37a) 
(375) 



Substitution into eq.(|T3|) yields the reduced generating functional 



1/2 



(x - X{0) 



+ 



1/2 



cosh^ 


[vl 


- 0(0) 


)1 




{x - X{0)) 


1 - ! (x - X{0)) 


(9 + Lf) 



(3^ 



which is still subject to minimization with respect to 0(o)(x), X{o){x) and L{x). Since 
the purpose of L is solely to impose the constraint / = 0, the above generating function 
is equivalent to: 



/ d^x'j 



1/2 



1 



{x - Xio) 



+ 



+ 



d^x'y 



1/2 



cosh^ 


.V 8 


(0 - 0(0)) 






{x - xm) 


"l - 


! (x - X(o)) 


9 



d^x7^/^ L{x) f (0(o)(a;),X(o)(a;) 
(minimized with respect to 0(o)(a;), X{o){^) L{x)) . (39) 

Here, I have set / = in eg . (|38D , and then appended a Lagrange multiplier term to it 
to recover the equivalent condition. 



6.1.1. Recovery of previous results The application of the Lagrange multiplier method 
to gravitational systems was suggested in an earlier paper by Salopek (1991). It will be 
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shown how these special results may be derived using the generalized methods presented 
in section 2. 

For two or more scalar fields denoted by 0a, a = 1, ■■■,n, a > 2, it was demonstrated 
in Salopek (1991) that one could obtain a special class of solutions to the long-wavelength 
problem by attempting an Ansatz which was an integral over the volume element of the 
Hubble function, H{(j)a): 

S[^abix),Mx)] = -2 / d'x^'/^H[Mx)] , (40a) 



provided that the Hubble function H^cpa) satisfies the separated Hamilton- Jacobi 
(SHJE): 

If one is fortunate to find a solution, 



H[<Pa\<Pa), (41) 

of the SHJE (^40^) which depends on n parameters, (pa, one may construct another 
solution H{(f)a), 

H{<f)a) = H{(pa\4>a), (42a) 

by choosing the parameters to minimize the n-parameter solution with respect (pa 
assuming that the parameters are subject to a constraint: 

/(0a) = 0. (426) 

By introducing a Lagrange multiplier A, one can write this new solution as 

H = H{<Pa\4>a)+Xf{4>a), (43a) 

where H is minimized with respect to A and (pa'- 

dH dH , 

Now in this paper, I have focussed on a single scalar field and a single dust field, 
but all past experience has show that a dust field can basically be treated as scalar field 
with a different term appearing in the energy constraint. In this subsection, for the 
purpose of simplicity, I will set g{(p(o), X(o)) = in eq.(^), and then obtain the following 
generating functional, 

4 f ,3 1/2 ^i^h' \\/li<P-<P{0) 



S=- fd'xj'/' 
3 J 



+ y rf^x7^/^ L(x) / (0(o)(x),X(o)(a;) 

(minimized with respect to 0(o)(x), X{o){x) and L{x)) . (44) 
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However, this generating functional is of the form of eq.( |40aD and eg. ( [43 oD where 



H = H 



0,xl0(o),X(o)] + A/ (0(o)(x),X{o)(a;)) , (A = -L/2) , 



is minimized with respect to A, 0(o) and X(o)- Here, H, 



H 



(p,x\(p(o),X(o) 



sinh^ 














X{0)] 





(45 a) 



(456) 



satisfies the following SHJE 



2 OH 2 fdH 

3 (9x 3 I 90 ^ 



(45 c) 



Hence the generalized Lagrange multiplier method recovers the special case considered 
in Salopek (1991) which was described for n scalar fields in eqs.( [43c| -b). 

I will now consider two special cases for / corresponding to initial hypersurfaces 
where firstly the dust field is uniform and where secondly the scalar field is uniform. 

6.2. Uniform dust field initially 

As was discussed before in section 2, if one takes 

/ (0(o)(a;), X(o)(a;)) = X(o) (46a) 

then variation of the Lagrange multiplier implies that 

X(o) = 0, (466) 

and that the initial hypersurface is one where the dust field is uniformly zero. In this 
case, I will further assume that g is an arbitrary function of • 



9 = 9i<P{0)) ■ 

One finds that the generating functional eq.(^) is: 



S[jab{x),(^{x)] 



tanh 



3 ^ X 



'X^ 



1 - zX9 



8 A, a0(oj 



(l - 1X9) ' 



(46 c) 

1 \ ,(47a) 
(476) 



The second equation is the minimization condition of the first with respect to 
For various choices of g, I have computed the resulting generating functional 
For g a constant: 



g = C (a constant) 



(48 a) 



For g proportional to an exponential: 

/3 



g = C exp 
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(486) 

(49 a) 
(496) 



For g proportional to cosh (^3/8 0(o)): 
g = C cosh"' 



4 /• , ,^ cosh^ (i/fi 



(50 a) 



(506) 



Thus the Green function method leads to many solutions of strongly coupled 
gravitational systems. 

6.3. Uniform scalar field initially 

Setting / = 0(0) implies that the scalar field is uniformly zero on the initial hypersurface. 
In this case, I now assume that g is an arbitrary function of X(o)- 

9 = 9 (X{0)) • (51) 
In this case, the generating functional eq.(|39D becomes: 

'^■^ [X- X(o) j 



tanh^ 



cosh^ (^y 






{x - X{0)) 




[i-t 




- xm) 


9 



x-xm)' 



9 



^ 16^ 



3 dg 



4 dxio) 



- I (x - X(o)) 9 

Equation (|52^) is the minimization condition of eg . (15201) with respect to X{o)- 
For various choices oi g, 1 have computed the resulting generating functional 
For g a constant: 

g = C (a constant) , 



(52a) 
(526) 



(53a) 
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5[7a.(x),x(x)] = C J dW^'exp 



(536) 



which is the same result found in eg . ([49 Bj) using a different route. 

For g = — 4/(3x(o)), one obtains 

4 1 . ^ 

5 = -o ' (54a) 



3 X(o) 

3 ^ ' X 



5 = _1 /'rf3^y/2l. (54^) 



For g = C/x{o), 



where 



C 

g = , (55a) 

X{o) 

A r , sinh^ + P) 

S = - Sx^'" ^^^^^^ . (556) 

3 ^ X 



sinh^/? = ^ . (55c) 
What is rather unusual is that for 

C = ~ (56a) 
one finds f3 = i7T/2, and the resulting generating functional is 



4 / .3^,1/2 (Vl 



d'x^'^' iX^, (566) 

3 ^ X 

which is certainly a solution to the energy and momentum constraints but it is not the 



solution given in eg. ( [54Q| ,b)! However, note that when = 0, they do agree — - they 
share the same initial functional, X. These two results reflect the multivalued solutions 
of the minimization eguation ( [52 6| ) . 



7. Lagrange multiplier method: gravity, dust, massless scalar field with 
cosmological constant 



I will repeat the analysis of section 6 but 1 will now add a cosmological constant to the 
system of gravity, dust and a massless scalar field. The Green function Q was given in 

eg-O- 
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When the functionals, X and £, 



(57a) 
(576) 



which determine the initial setting do not contain spatial gradients, then the 
minimization of S, eq.(O), with respect to the initial 3-metric is straightforward: 



giving 



z = 0. 
1^7(0)' 

S 

Wo) 



1/4 



cosh 








cosh 


[Ma 


[x - X(0)) 




^'^hJ smh 







(58 a) 
(58&) 



/ o 1/2 r3-ffo / \ 
d a;7 / cotanh ^- (^x - X{o) j 



+ 



d^X'j 



{2Ho 



1/2 



cosh 



sinh 



^ (X - X(0)) 



cosh 



^ (X - X(0)) 



2^^ sinh 



M^ 

2 



d'x^'^^f. 



^ (x-X(o)) 
(59) 



In the above equation, I followed the analysis of section 6 by setting / = in the 
generating functional and then recovering this same condition by adding a Lagrange 
multiplier term which is a linear functional in L{x). 

7.1. Uniform dust field initially 

I will assume that initial state is given on an initial hypersurface where the dust field is 
zero which is characterized by a vanishing value of X{o){^)'- 

f = X(o){x) = 0. (60) 

The Lagrange multiplier term will thus be removed. Moreover, I will assume that g, 
which characterizes the initial state, is a function soley of the initial scalar field 0(0): 

9 = 9 (0(0)) • (61) 
After minimizing with respect to 0(o), the generating functional is: 

5[7afe(x),0(x),x(x)]/(2ifo) = - Jd^x-f^/^ cotanh^ + 
1 cotanh^ — h 



I 



d'^X'j 



1/2 



sinh 6 r 



cotanh^^ — h (0(o)) 



2 f dh 

3 \^ci(/>(o) 



2 ' 



(62 a) 
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where 9 and h are dimensionless representations of x ci'iid g, 



9 



h = g (0(0)) /{2Ho) . 



(626) 



The parameter field (p{o){x) is determined imphcitly through the algebraic equation, 

^ dh _ 

(62 c) 



tanh 



3 (i(/>(o) 



cotanh6' — h 



7.1.1. Exact solution One obtains a non-trivial exact solution if one assumes that g 
has the following form: 

9 (0(0)' 



2Hn 



= h = Cexp\ -y^^io) 



D + Eexp I \i-(P[o) 



(63 a) 



One can then solve 0(o) explicitly using eg. ( |62c| ), and the generating functional is 

S[^ab{x)A{x),x{x)] = -2 J d'x^'^^H{(f),x), (636) 
where the Hubble function H{(j), x) is given by 
H{(f),x)/Ho = cotanh^^ 

n -I- 

(63 c) 



1 [Eb + C/b - D + cotanh^] 



sinh'^ [(cotanh^ -D)'^-ACE] 



and b denotes 



exp 



As a check of this method, one may verify that if is a solution of the SHJE, 



H' = -- 



2 OH 2 fdH' 



+ - 



-I 



3 5x ■ 3 ^ _ 

Some special cases illuminate the solution. 
li D = C = E = 0, the Hubble function is 

f3HoX\ 



(63 d) 



(64) 



H = Hq tanh 



V 2 / 



, (^D = C = E = 0). 



li D = 0, C = E 



-|, the Hubble function is 



H = Hq cosh 



[D = 0, C = E = -- 



(65) 



(66) 



li D = and C = E, the Hubble function is 

H (1 - 4:E^) cosh^ sinh^ - 2E cosh , ^ ^ 



cosh^^ - 4^2 sinh^^ 



{D = 0, C = E) . (67) 
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8. Advanced example: initial hypersurface condition containing spatial 
gradients 

I will now consider an advanced example where the Lagrange multiplier term contains 
spatial gradients. This example is particularly important because it illustrates how the 
energy constraint applies to the initial fields. 

The system under examination will contain only gravity and dust, and the Green 
function is given by eq.(^^c). The function / which defines the initial hypersurface will 
be chosen to be, 



/ = X(o) - BR, 



(0) , 



(68 a) 



where -B is a constant and -R(o) is the Ricci scalar associated with the initial 3-metric. 
The functional X is taken to be the simplest non-trivial example, namely the volume of 
the initial 3-geometry: 



X = A I d"x 7(Q) 
The full generating functional is the sum of Q, X and C: 

[2y/Si^;cosh(./f.)-7^/^-7;o) 



(685) 



- / d^x 
3 



1 



3 X 1/9 1/2 



1(0) H^) 



-A I d^x-fli^+ I d^x-^'^ 



X(o)(a;) - BR(o){x] 



(minimized with respect to 7^°''(x), X(,o){x) and L{x 



(69) 



This example will not be exactly solvable unless B vanishes. One obtains approximate 
results by expanding in powers of B. 



8.1. Minimization conditions 

Variation with respect to L{x) and X(o) lead immediately to: 

= X(o) - -Bi?(o) , 



L = - 



4 



0)) cosh (yfz) - (7/7(0))" 



(x - X(o) 



(70 a) 
(70 &) 



After minimizing eq.(^) with respect to 7^^^ the trace and traceless parts are (the trace 
is defined using the initial 3-metric): 



= 2 



(7/7(0))^^^ cosh(yf2; 
X - X(o) 



Q /I 

— + 5 (i?(o) L + 2L^%) , (70c) 
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= 6-^/2 



1/4 



sinh 



X - Xio) 



In ([h-'] [/.(o)] 



{70d) 



8.2. Initial energy constraint 

After one minimizes 5, the initial fields, 7^°^(x), X(o)(a^) and L{x), are related in a subtle 
way. By construction, the Green function Q satisfies the energy constraint: 







+ 7 (27„c Ibd - lab led) 



60 50 



(71) 



()x(a;) (37ab(x) ()7cd(x) 

We interpret this equation as stating that the original fields satisfy the energy constraint. 
However, because of a symmetry transformation that relates the original fields and the 
initial fields, 



{lab{x), 



X{x, 



(7l?(x), 



-X{o)ia;, 



the initial fields, lfi^{x) and X(a){'x)i also satisfy the energy constraint: 







50 



'^(0) \^lac Ibd lab led 



60 



50 



^labi^) ^lZ\x 



(0), 



(72) 



(73) 



SX{o)ix) 

By virtue of the minimization conditions (see,e.5'., eq.(^6^d)), one may thus replace 
functional derivatives of with functional derivatives of (X + £) , 

60 S(J + C) 



6X{o){x) 6X(Q){x) 
60 6(1 + C) 



57l?(x) 57l?(x) ' 
in eq.([73|) which leads to the initial energy constraint, 



(74a) 
(746) 







6{I + C) 



^(0) 



1/2 



'^lac Ibd lab led 



5(J+£) S{I + C) 



(0), 



(75) 



6X{o){x) 
(initial energy constraint) , 

which is independent of the original fields, iab{x) and xi^)- Explicit evaluation of the 
functional derivatives in eq.(|75D yields a relationship amongst the Lagrange multiplier 
and the initial fields: 

Q = L + 2B'^rrf^mab--^ " ^ 

3 

where the tensor m""^ is given by 



(initial energy constraint) (76a) 



ab 



Rio) L + 7(0^) L 



ab T ;c 



:ab 



(766) 
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m 



ab 



with m 



^(0) ^ab 
lab '''' 



(76 c) 



and denotes its traceless component: 

1 

3' 

In the above, a semicolon (; ) denotes a derivative with respect to the initial 3-metric 
and all indices are raised and lowered using the initial 3-metric. In principle, one may 
solve for the Lagrange multiplier in terms of the initial 3-metric, although this is difficult 
to achieve in practice since this is a fourth order, nonlinear partial differential equation. 



.3. Classical evolution 



Classical evolution is found by solving eqs.( [70a| -d) for 'jab{x) in terms of the proper time 
x{x), the initial fields, 7*^°^(x) and X{o){^)^ a^iid the Lagrange multiplier, L{x): 



7 
.7(0), 



1/2 



[h] = exp 2z 



In 



{x - X(o)) 

X(i) - X{0)] 

-X(0)) 



(x - X(o)) 

X(2) - X(0) 



1 



[X(2)-X(0)J 

(X-X(O)) 
[X(l)-X(0)) 



(77a) 
(776) 

(77c) 



^ is a quadratic function of (x~X(o))- The terms, — X{o)) and 

X(2) - X(o)j, denote the two roots of (^7/7(0) j : 



(7/7(0))" 



X(i) - X(o) 



A 25 / 2 ^.^ 

A IB 



L+ 3 



/?(0) + - V'-.c ) + 



/85 
3L 

3L 



{lid) 



{lie) 



,X(2) -X{0), 

The tensor nf"^ was given earlier but I will write it down here just to have all of the 
results in one place: 

m'^" = R$-^ L + 7^') - i^'"' . (78) 

X(o) is given by the initial hypersurface condition: 

= X(o) - 5i?(o) , (79) 

and the Lagrange multiplier may be solved in terms of the initial 3-metric from the 
initial energy constraint: 







L + 2B'^rrf^ m^fe - ^ (^5 m + '^A 



(80) 
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One may solve for L, and subsequently all other quantities, by expanding in powers 
of S. To first order in B, one finds: 



(l + BAR^o)) , 



X(o) = BR(^o) , 
^3A^\ 



m 



m 



ah 



V 4 

'3^2^ 



-R(o) + BA (i?(o) i?('o) - -R(o) + 7(*o) D'^R{o)) 
R{o} + -B^ (-^(0) + 2-D^-R(o) 



4 \ / 

— 111 :— -ft{0) 



X{i) - X(o) = ^- 2 



-Rem ) - W - -B ( Rlo^Rab"^ 



X(2) - X(0) 
35^2 



V3a; 

X 



Rio)]+\I^B [rIr^'^ 



1/2 



1/2 



-pa6 -p(0) 

2 /•. 3Ay\ \^{0)^ah 



lab 



(0) 
Tab 



1/2 



(81a) 
(816) 

(81c) 
(81 (i) 
(81 e) 

(81/) 
(81^) 



(3^xi?i? + (l-^)i?(o)7, 



In the above, the laplacian of i?(o) is denoted by 
D^R{o) = (-R(o))' ;c- 

<5.^. Determination of the classical 4-''netric 



/aft 



(82) 



In solving the classical Einstein equations, one ordinarily computes the 4-metric 
describing time and space. However, in the HJ formalism, one's attention is primarily 
focusscd on the 3-mctric which describes the spatial geometry. How docs one recover 
the 4-metric? In general, given the generating functional S['yab{x),x{^)]j 4-metric is 
computed by making an arbitrary choice for the time parameter, and then integrating 
the definition of the momenta, 

5S 



(x - N^X,i) IN = K, 



(83 a) 
(836) 



(81/i) 



which is valid in the strongly coupled limit. 
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For the problems considered in this section, a natural choice for the time 
hypersurface is one where the dust field x is uniform, 

t = xl^ (84a) 



which describes a comoving slice. If one assumes that the shift Ni vanishes, then eq.(^3^) 
implies that the lapse N is equal to one, 

AT = 1, (84&) 

which describes a synchronous gauge. The line-element describing the 4-geometry is 
then 

ds^ = — ^ + 7afe rfx" dJ' , (85) 

where the classical evolution of the 3-metric ^ah was computed in the previous section 
using the HJ minimization prescription; see eq.( |81/j| ) for the first few terms. At each 
spatial point, the 4-geometry is locally Kasner (see, for example, Salopek 1998). 

8.5. Computation oj generating junctional 

Using the exact expressions for classical evolution, eqs.( |77o| -e), one may express the 
generating functional S, eq. (|69D , in an elegant form which depends on the dust field 
x(a;), the initial fields and the Lagrange multiplier: 

S = -\ d'x^li^ l{x- Xm) + X + £ . (86) 

Here I have eliminated all reference to the original 3-metric '~fab{x). As was found in 
Salopek (1998), the generating functional is basically linear in x(x). If L{x) is positive, 
the integrand of S decreases in x(x), whereas the opposite is true if L{x) is negative. 

Ultimately, one wishes to compute S soley as a functional of the original fields, 
labix) and x{x)- Since this cannot be done exactly for the problem at hand, I will be 
content to approximate it using a Taylor series in B. First note that the initial 3 metric, 
7^°''(x) may be expressed as a function of "jabix) and xi^)'- 



lab — '^ab 



(87a) 



which is accurate to first order in B. Here the tensor kab is conformally related to the 
original 3-metric, 7afe, 

fcafc= (l-^) "\ab. {87b) 

and R'', i?^^ denote the corresponding Ricci scalar and tensor, respectively. After 
imposing the initial hypersurface condition, eq. ([75|) , one may safely drop the Lagrange 
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multiplier term L in eq.([86|) and the generating functional becomes, 



4 



which is accurate to first order in B. Conceptually, there is no problem in extending 
this expression to higher order in B. 

9. Summary and conclusions 

A strong coupling expansion appears in many investigations of gravitational systems 
including long-wavelength cosmological fluctuations, gravitational collapse and string 
theory formulations of cosmology. In the present paper, powerful Hamilton- Jacobi 
methods have been further developed to allow for a very general solution of strongly 
coupled gravitational systems. 

In section 2, one generalizes the semiclassical Green function method of solving the 
constraint equations to encompass the specification of an arbitrary initial hypersurface. 
One assumes an Ansatz for the generating functional 5, eq. (p!3D , which is the sum of 
three terms: 

S = g+X + C. (89a) 
The Lagrange multiplier term. 



C = j d^x^l^^ L{x) f [^ab'i^)^ <Pio){x), X{o){x)\ , (896) 

specifies the initial hypersurface, 

/ = 0. (89c) 

I is the initial state. The Green function Q determines how the systems evolves from 
the initial setting. The Ansatz (|89q|) is justified mathematically by computing its 
functional derivatives. It hence satisfies the energy constraint. Because gauge-invariance 
is maintained at each step, the Ansatz satisfies the momentum constraint. Classical 
evolution follows from minimization of the generating functional ( ^9q| ) with respect to 
the initial fields and the Lagrange multiplier L. 

To illustrate the generalized method, I constructed in sections 4 and 5 Green function 
solutions for (1) gravity, dust, a massless scalar field and a cosmological constant and 
(2) gravity interacting with a scalar field with exponential potential. One verifies these 
solutions after deriving a reduced energy constraint in section 3. 

If the functionals, X and C, which define the initial setting do not contain spatial 
gradients of the initial fields, then the program is straightforward to implement, and one 
recovers a primitive form of the minimization principle that had been advanced in an 
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earlier paper (Salopek 1991). In fact, in sections 6 and 7, exact solutions were derived 
that had not been previously known. 

However, in general the initial hypersurface condition ( |89 c| ) contains spatial 
gradients, and one must typically resort to an approximation method as was illustrated 
in section 8 for the case of gravity and dust. There it was shown explicitly how to deal 
with the initial energy constraint. 

The semidassical Green function method has reached a high level of generality which 
should be sufficiently powerful to treat most strongly coupled gravitational systems of 
physical and numerical interest. Hamilton- Jacobi methods have been useful in the past 
for solving nonlinear problems in cosmology. In the future, the methods presented in 
this paper may shed some light on quantum aspects of the gravitational field. 
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